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A central problem in gravitational wave research is the generation problem, i.e., 
the problem of relating the outgoing gravitational wave field to the structure and 
motion of the material source. This problem has become, in recent years, of in- 
creased interest in view of the development of a worldwide network of gravitational 
wave detectors. We review recent progress in analytical methods of tackling the 
gravitational wave generation problem. In particular, we describe recent work in 
an approach which consists of matching a post-Newtonian expansion of the met- 
ric near the material source with a multipolar-post-Minkowskian expansion of the 
external metric. The results of such analytical methods are important notably 
for providing accurate theoretical predictions for the most promising targets of 
the LIGO/VIRGO interferometric network: the "chirp" gravitational waveforms 
emitted during the radiation-reaction-driven inspiral of binary systems of compact 
objects (neutron stars or black holes). 



1 Introduction 

I wish to dedicate this talk to Henri Poincare who introduced (several aspects 
of) the concept of gravitational wave ("onde gravifiqne") ninety years ago. 
Indeed, in his two seminal papers of JuneB and JulyEl 1905 (the first one of 
which preceded Einstein's paper on special relativity by one week), Poincare 
not only introduced what was to become later the basic mathematical struc- 
tures of special relativity (the Poincare group and the "Minkowski" metric 
{it) 2 + x 2 + y 2 + z 2 ), but also pioneered the idea that one needed, for consis- 
tency, a relativistic theory of gravitation. In his two papers written in 1905, 
he defines a class of "Poincare invariant" gravity theories and emphasizes that 
they .predict that the gravitational interaction propagates with the velocity of 
lightu. In a later work of 1908a, he went as far as speaking of the emission of 
gravitational waves ("onde d'acceleration") and of the associated loss of energy 
of the emitting system. He even mentions that the main observable effect of 
this dissipation of energy into gravitational waves will be a secular acceleration 
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of the mean motions of planetsu. It is interesting to note that .the discovery 
of binary pulsarsB, and their subsequent continuous observation 0, has allowed 
one precisely to verify the aspects of gravitational waves discussed by Poincare: 
propagation of the gravitationaLinteraction with finite velocity, and associated 
effect on the orbital dynamics El. This verification has been observationally 
possible in two binary pulsar systems: PSR1913+160 and PSR1534+12§. 

The discovery of binary pulsars is important in three respects for gravita- 
tional radiation research: (i) it establishes the reality of gravitational radiation 
by verifying in a direct manner that gravity propagates with the velocity of 
light between the companion and-th^ pulsar; (ii) it gives us our first tests of 
the strong-field regime of gravity E2I, EU thereby confirming the validity of Ein- 
stein's theory in a regime so forth untested; (iii) it establishes the existence of 
strong sources of gravitational waves, thereby providing fascinating targets for 
the LIGO/ VIRGO network of interferometric detectors. 

Indeed, the observation of the secular acceleration of the orbital mean 
motions of the binary pulsars 1913+16 and 1534+12 proves that, in a few 
hundred million years, these systems will have shrunk so much that they will 
constitute an "inspiralling" binary system of neutron stars: i.e. a very close 
system of two neutron stars, orbiting around each other at a very fast and 
accelerated pace, the orbital frequency increasing from, say, ~ 10 Hz to ~ 1000 
Hz in about 20 minutes. Then, when the two stars get near each other they 
start coalescing together to form only one central object. In the last minutes 
of the inspiralling motion such systems emit rather strong gravitational waves. 
The first person to conceive of such systems, and to realize they provided 
superb targets for gravitational wave detectors, was Dyson E3, in a prescient 
paper written years before the discovery of pulsars established the existence 
of neutron stars. The characteristics of inspiralling and coalescing binari£s_as 
graiatational wave sources have then been explored by several authors Ej,li-3, 
Ejjlij. The rate of occurrence of such events seem to be high enough to fuxnialr 
a regular (~ monthly) source of signals for the LIGO/VIRGO network E3, Ea. 
The information content of inspiralling events is of excellent quality: their 
detection (if done with a suitably high signal to noise ratio) should allow one: 
(i) to measure directly cosmological distances and. i^r^y_to_have a clean 
access to the cosmological parameters H and q Ej, Ej, E3, Ej, l3; (ii) to test 
the nonlinear structure of radiative gravity E3; pii) to perform new tests of 
the existence of a scalar component to gravity Ej; (iv) to probe black hole 
physics Ej. In view of the importance of inspiralling events, it is crucial to 
dispose of an accurate theoretical model of the corresponding gravitational 
wave signals. The aim of the present contribution is to review recent progress 
in analytical methods of tackling the generation of gravitational waves, and 
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their application to inspir ailing binaries. 

2 Analytical formalisms for treating the generation of gravitational 

waves 

One can distinguish three basic questions in gravitational radiation theory: 

- Question 1 ("asymptotic problem"): What is the asymptotic behaviour, 
appropriate to isolated systems and consistent with Einstein's field equations, 
of radiative gravitational fields far from their sources? 

- Question 2 ("generation problem"): What is the link between the preced- 
ing asymptotic behaviour and the structure and motion of the sources that 
generate the gravitational radiation? 

- Question 3 ( "radiation reaction problem" ) : What is the back-reaction of the 
emission of gravitational radiation on the source? 

The standard answer to Question 1 is given by the Bondi-Sachs-Penrose de- 
scription of radiative, asymptotically flat spacetimes, with a sufficiently smooth 
fall off at 2 + and the exclusion of ingoing waves on I~ . However, this answer 
is still a conjecture. Let us note, in this respect, that^ the estimates used in 
the global theorem of Christodoulou and Klainermantj are not strong enough 
to establish the standardly assumed peeling at I + ; and (ii) some perturbation 
calculations suggest a violation of peeling in scattering problems E3 (where it 
is found that I + cannot be C 3 ). 

Questions 2 and 3 have standard answers (discussed in textbooks) only 
at the lowest approximation. These standard answers go by the (ambiguous) 
name of "quadrupolc formulas" . Actually, one should carefully distinguish: 
the "far-field quadrupole formula" , the "energy-loss quadpipole formula" , the 
"radiation-reaction quadrupolar force", etc... (see, e.g.,c3 for a discussion). 
In any case, these standard answers are insufficiently accurate to give mathe- 
matical models of inspiralling signals adequate for high-precision observations. 
Indeed, during the final stages of the inspiralling motion the orbital veloci- 
ties become rather high (u/c<0.3) and necessitate the consideration of many 
corrections to the leading "quadrupole" result. 

Several methods have been proposed for going beyond the loweat-jHrdcr 
results. For instance, some years ago Epstein, Wagoner and ThorneEj, Ej (in 
an attempt to generalize the Landau-Lifshitz-type derivation of the standard 
far-field quadrupole formula) introduced a post-Newtonian extension of the 
quadrupole formalism. Though their formalism is marred by some mathemat- 
ical difficulties (divergent integrals), it was used to-derive 0(v 2 /c 2 ) correc- 
tionsL-tp the quadrupole formula for binary systems Ed, and has been recently 
usedE3 to go to higher orders in v/c. A basic problem of the Landau-Lifshitz- 
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Epstein- Wagoncr-Thornc approach is the lack of a clear separation between 
the near zone and the wave zone. The combined use of an "effective" stress- 
energy tensor for the gravitational field (with non-compact support) and of 
formal post-Newtonian expansions quickly leads to the acmeamnije jaf .diver.- 
geiiL integrals. By contrast, Blanchet, Damour and Iyer Ej, O, £3, E3, E3, L3, 
cB, c3 (building on the Fock-type derivation, of the quadrupole formula and on 
the double-expansion method of Bonnor EH) introduced a new gravitational- 
wave-generation formalism based on a clean separation between near-zone 
and wave-zone effects. The Blanchet-Damour-Iyer approach is mathematically 
well-defined and obtains corrections to the leading quadrupolar formalism in 
the form of compact-support integrals. [In a recent development of this for- 
malism, BlanchetEa found it convenient to obtain the corrections in the form of 
(well-defined) analytically-continued integrals which are (formally) equivalent 
to compact-support integrals.] The BDI scheme has a "modular structure" : the 
final results are obtained by combining an "external zone module" (in which 
the external, vacuum metric is expanded as a multipolar-post-Minkowskian 
double series) with a "near zone module" (based on a more traditional post- 
Newtonian-type expansion) . When dealing with strongly self-gravitating mate- 
rial sources (such as neutron stars, or black holes) one must also use a "compact 
body module" c3. After elimination of the various mathematical intermediaries 
appearing in the formalism (such as the "algorithmic" multipole moments Ml 
and Sl), the basic structure of the final results of the BDI formalism is the 
following: the (directly observable) "radiative" multipole moments Ul and 
Vl Cj, parametrizing the angular dependence of the asymptotic gravitational 
wave amplitude /i^ T (T, R, 6, ip), are given in terms of the "source" multipole 
moments II and Jl as a series of terms of increasing nonlincarity: 

U L (t)=d*I L (t)/dt*+F L lI(t'),J(t% 

v L {t) = d e ,j L (t)/dt e + g L [i(t'), j(t% 

where Tl and Ql are multilinear (successively quadratic, cubic, etc. . .) retarded 
functional of the full past behaviour of the source moments Il> it') , Jl 1 it') (for 
f < t). As for the "source" moments I Lit), J Lit) they are (approximately) 
defined in the formalism as some instantaneous (nonlinear) functionals of the 
source variables ^"(x, t) [or rriA, in the case of spherical compact bod- 

ies!. 



For instance, the first two corrections.-(p 
radiative quadrupole moment are given byL3, 
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with 



h(t) = J<Pxx L a(t, x) + m I 3)c2 f*±S***Q, x) 



4(2^+1) 



{£+ 1){2£ + 3)c 2 dt 



d 3 xx kL a k (t,x) + 



(2) 



where a = (T°° + T kk )/c 2 , a % = T M /c. The integral appearing on the right 
of equation [j] represents the effect of the backscattering of the gravitational 
waves on the Sckwaczschild-like curvature associated to the total mass M of the 
source ( "tails" ) E-a, E3. Beyond the terms written in equation [l] there are many 
other nonlinear contributions (of formal higher order in v/c). For instance at 
the quadratically nonlinear order one has a contribution to f/r,(i) depending 
upon the gravitational wave flux emitted in the past, 

^-H\ f f dE^ w (t') \ 
(£+l)(l + 2)J_J t { dt'dn ) L ' [6) 

which has beemilisims^dJin. different guises, and under different names) by 
several authors ej, ^q3-Cj-q Fiir^plicit. applications of the presently dis- 
cussed scheme, see EJ, Ell, Ej, Ej, o, O, tJ, lj, EJ. Besides "hereditary" effects 
in the wave zone (such as the integral contributions in equations |l| and |^), it 
has also been possible to investigate the leading hereditary effects appearing in 
the near-zone field: they enter at the fourth post-Newtonian level, i.e. (v/c) 8 
beyond the Newtonian approximation, and correspond to "rtail" modifications 
of the 0(v 5 /c 5 ) Burke-Thorne radiation reaction potential E3. 



3 Inspiralling compact binaries 

The accurate mathematical modelling of the gravitational wave signals emitted 
by inspiralling compact binaries needs two (related) inputs: (i) a solution of the 
"generation problem" ; and (ii) a solution of the "radiation-reaction problem" 
adequate for treating compact objects. The schemes discussed in the previous 
section were primarily aimed at solving the generation problem, i.e. at giv- 
ing ft,^ T as a retarded functional of the structure and motion of the source. 
However, to have an explicit representation of hf^ as a function of time, one 
needs to know the time evolution of the source, i.e. to solve the problem of 
motion, including radiation-reaction effects. Actually, it has been recently em- 
phasizedc3 that the radiation-reaction part of the problem was the most crucial 
one in that it determined the time evolution of the phase of the gravitational 
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wave signal (which follows, modulo a factor two for circular orbits, the orbital 
phase). [We work here within the "restricted waveform" approximation, i.e. we 
focus on the main Fourier component of the signal.] Indeed, an accurate mod- 
elling of the (radiation-reaction-driven) phase 4>gw (*) of the gravitational wave 
is essential for a successful detection based on correlating the observed, poiH4 
h obs (t) with some theoretical template h thcmy (t) = a GW (t) cos <f>Gw(t) 0,0,0 
(more so than the modelling of the evolution of the amplitude aawiff)- The 
results discussed in the previous section are certainly accurate enough for pre- 
dicting cicw (t) ■ The situation for what concerns the radiation-reaction driven 
phase <t>Gw{t) is less satisfactory. Indeed, the only complete results available 
for the equations of motion of a compact binary are the (v/c) 5 -accurate equa- 
tions of motionB. Beyond this levelj-anlv partial results are known: namely the 
0((v/c) 7 ) radiation-reaction terrpjE2l,E2l, and the 0((v/c) 8 ) hereditary contri- 
bution to the radiation-reaction E3. The (expected) "conservative" contribu- 
tions (of order (v/c) e + (v/c) s + • • •) to the equations of motion are unknown, 
as well as the higher-order contributions to radiation reaction. The only way 
one can presently deal with this problem is to heuristically rely on a (naive) 
energy-balance argument to relate the loss of mechanical energy of the binary 
system to the asymptotic flux of gravitational waves. In technical terms, one 
writes 

^gmcchanical dE^^ ^ UX 

= X/ r 2£+l 



dt dt 

e=2 



(£ + !){£ + 2) f dU L \ 2 41(1 + 2) f dVL \ 

(e - i)m(2e + iy.\ \ dt J + c 2 (i-i)(£ + iy.(2£ + iy\ \ dt J 

(4) 

in the right-hand side of which one inserts the best available results (from 
solving the generation problem) on the radiative multipole moments generated 
by a compact binary. 

The highest-accuracy results obtained along these lines have been derived 
in the limiting case of a very small test mass orbiting a heavy central mass 
modeled as a Schwarzschild (or Kerr) black hole. Indeed, in such a case one 
can treat the generation problem as a linear perturbation of Schwarzschild or 
Kerr. The latter pertubation problem can, thanks to the work of many people 
(Regge, Wheeler, Zerilli, . . ., Teukolsky, . . ., Sasaki, Nakamura, . . ., Chan- 
drasekhar, ...), be reduced to integrating some linear ordinary differential 
equations, for the radial dependence. This integration can be done numeri- 
callyEa, E3. Moreover, the post-Newtonian expansion (in powers of v/c) of the 



G 



generated gravitational wave amplitudes can be derived analytically S, El, El. 
The results of this approach are important testbeds for the full problem and 
can help us in posing the important questions (such as: "how fast does the 
post-Newtonian expansion converge?" E3) , but they fall short of providing us 
with the answers we really care about. Indeed, if we introduce the dimension- 
less parameter measuring the deviation from the test mass limit, 

mim 2 ... 

v = ^ = 7 — ; — ^ ( 5 ) 

(mi + m2) 

(where m% and mi are the masses of the members of an inspiralling binary), we 
expect that many of the systems that LIGO /VIRGO will detect will be made 
of two nearly equal neutron star masses: m\ 7712 ~ 1.4M Q . This corresponds 
to the maximum possible deviation from the test-mass limit [y = j instead of 
v <C 1) for which the black-hole perturbation results become unreliable. This 
is why we badly need the general-purpose analytical methods discussed in the 
previous section. Only such methods can, at present, deal with inspiralling 
binaries having comparable masses m\ ~ mi- [Note, however, that black-hole 
perturbation methods are directly relevant for dealing with some of the target 
sources of low-frequency space interferometers such as LISA: e.g. the fall of a 
neutron star into a very massive black hole.] 

The main result one is interested in is an analytical expression giving the 
time evolution of the gravitational wave phase, i.e. something we can call the 
"phasing formula" of inspiralling binaries: 

(j) GW = ^ORBITAL = F[te . pzl (g) 

where t® is the (proper) time at the Earth laboratory recording the continuous 
arrival of the (main) gravitational wave signal 

h% w {t (B ) = a GW (t e ) cos ^ GW (t e ), (7) 

and where {pi} is a set of parameters carrying information about the emitting 
binary system. It is interesting to note that the "phasing formula" is nothing 
but a continuous analog of the discrete "timing formula" which is bai 
relativistic pulsar timing. Indeed, the timing formula of binary pulsars 
can be written as 

4>% SR =F[t®; Pi }, (8) 

where t% is the (proper) time of arrival at the Earth observatory of the iVth 
(N G N) pulse emitted when the rotational phase of the spinning pulsar was 
(f)PSR ^ 27riV+const. Here also {pi} is a set of parameters carrying information 
about the binary system, and the basic aim of pulsar timing is to measure as 
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many p^s as possible 0,0. The analogy between equations |6| and ^ is clear: 
in one case one is timing from Earth a continuous orbital phase, in the other 
case one is timing a stroboscopic rotational phase. 

Thanks to recent works which explicitly cWekinfid the analytical methods 
discussed in the previous section to the (v /c) 4 |£3,E3,lj and (w/c) 5 c2l orders one 
knows the phasing formula of inspiralling binaries (of arbitrary masses) to the 
following accuracy: introducing the dimensionless time variable 



(?V 



5G(toi + m 2 ) 



(where t c denotes the coalescence time) we can write 

1-5 



iORBITAL 



-ft|l + A 2 (i/)T~i + Az{v)t^i + Ai{v)t~i 

, N . (10) 
■B 5 (v)t-slnt + 0[t-i) |, 



where 

. , , 5 /743 11 , 

MW) = — 1 — v (n) 

■ ' 24 V 336 4 ' 



(v 2 /c 2 corrections: see0,E3) 



Mv) = -\* (12) 



(v 3 /c 3 or "tail" corrections seeEl for the v = limit, and0,0,El for v ^ 0) 

5 / 1855099 56975 371 



v; 64 V 225792 4032 32 / v ; 

(v 4 /c 4 corrections: seeE3 for the v — limit and§,ll,@ for v 7^ 0), and 

. . / 38645 15 \ 

Bs{v) = -it (^^72032 + 2048 / 

(xr lj£ corrections, after the absorption of any A 5 (v) in the definition of 4> c : 
see l£j for the v = limit and u£l for v ^ 0) . The numerical importance of 
finite-mass-ratio effects (y ^ 0) is to be noted. In particular, the equal mass 
case (y = jj represents (with respect to the test-mass limit) an increase of Ai 
by 31% and of A4 by 52%! (By contrast the corresponding change in B§ is 
only 0.8%.) 
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4 Conclusions 



• A general comment on the present brief review of analytical approaches to 
gravitational radiation is that it confirms the perennial validity of a remark by 
Poincare: to the effect that real problems are never definitively solved, but only 
more or less solved ("il y a seulement des problemes plus ou moins resolusEa). 

• Though I have insisted on the need of computing the (hard to get) higher- 
order contributions to the phasing formula of inspiralling binaries because of 
their importance for extracting the maximum possible information from gravi- 
tational wave signals, it should be stated that these corrections are (probably) 
not needed for searching and discovering gravitational wave signals in the noise. 

• The ultimate post-Newtonian accuracy which is really needed in the phasing 
formula for an acceptably accurate determination of the information-^txypg 
jPjirameters {p{\ — {t c , v^(m\ + 7712), v} is still unclear at present (seeE3,E2l,Ea, 

• I have considered above only the simplest case where the members of an 
inspiralling system are slowly spinning. This is the .situation one can plausibly 
expect in most neutron star- neutron star systems E3. However, systems con- 
taining black holes (if theji exist in appreciable number) might contain fast 
spinning objects. See e.g.E3,EZI for spin-dependent effects. 

• If very high post-Newtonian contributions to the phasing formula arc really 
needed, one might need to reconsider the presently developed analytical ap- 
proaches. It might, for instance, become necessary to define them in a fully 
algorithmic manner allowing the use of computer-based algebraic programmes, 
or it_might become necessary to match them to numerical relativity results (see 
e.g. lid, L5). A better understanding of the mathematical nature of the post- 
Newtonian expansion might also help. 

• I anticipate that a serious obstacle to improving the present (v/c) 5 accuracy 
of the phasing formula will come from the need to extend the accuracy xif the 
equations of motion of compact binaries beyond the G 3 level treated inll3. 

• The increasingly slow convergence of the post-Newtonian series toward the 
end of the inspiralling stage points out the importance of improving the sensi- 
tivity of the detectors on the low- frequency side (say a few tens of Hz, corre- 
sponding to gravitational waves emitted early on). I note, in this respect, that 
the VIRGO detector puts a particular emphasis on improving its low-frequency 
sensitivity. 
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